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Abstract 

The recently formulated framework of anisotropic hydrodynamics is used in 3+1 dimensions 
to study behavior of matter created in relativistic heavy-ion collisions. The model predictions 
for various hadronic observables show that the effects of the initial anisotropy of pressure may 
be compensated by appropriate adjustment of the initial energy density. In this way, the final 
hadronic observables become insensitive to the early stage dynamics and the early thermaliza- 
tion/isotropization puzzle may be circumvented. 



1. Introduction 

Soft-hadronic observables measured in the relativistic heavy-ion experiments are well repro- 
duced by perfect-fluid hydrodynamics or by viscous hydrodynamics with a small viscosity to 
entropy ratio |lll|2|,|3l|4|,|5l|6|]. Nevertheless, the use of such approaches at the very early stages of 
the collisions encounters conceptual difficulties. Thermalization times shorter than a fraction of 
a fermi (used in the perfect-fluid approaches) cannot be explained within microscopic models of 
the collisions. On the other hand, viscous hydrodynamics is based on an implicit assumption that 
one can make an expansion around an isotropic background. If the shear correction is large, a new 
framework incorporating large momentum-space anisotropics into the leadingorder of the ap- 
proximation may be useful. Such a new approach has been introduced in 10, SM [13; 11. 12. 131 
and we refer to it below as to the anisotropic hydrodynamics. 



2. Anisotropic hydrodynamics 

Anisotropic hydrodynamics is based on the equations 

df.T^"' = 0, (1) 
d^cT^ = £, (2) 

which express the energy-momentum conservation and entropy production laws. The energy- 
momentum tensor T^^ has the structure 

T^"' = (e -H U^'U" - - (P ^ - P\\)WY , (3) 
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where f n is the longitudinal pressure and is the transverse pressure. In the limit P\\ - P± = P, 
Eq. (O reproduces the energy-momentum tensor of the perfect fluid. Similarly, the entropy 
production law (|2]i is reduced to the entropy conservation law, if we take I, -Q. The four-vector 

describes the flow of matter, while V defines the beam (z) axis. In the general case, we use 
the parameterizations = (mq cosh §, u^, «,„ mq sinh §) and = (sinh §, 0, 0, cosh)?), where 
and My are the transverse components of the four-velocity field and § is the longitudinal fluid 
rapidity. The entropy flux cr'^ equals cr U^, where cr is the non-equilibrium entropy density. 

One can show [JJ that instead of Py and P± it is more convenient to use the entropy density 
cr and the anisotropy parameter x as two independent variables (one may use the approximation 
P\\/Pi. ~ Jc"^^"*)- Similarly to standard hydrodynamics with vanishing baryon chemical potential, 
the energy density e, the entropy density cr, and the anisotropy parameter x are related through 
the generalized equation of state |@] 

s(x,o-) = eqgp(cr)r(.x), (4) 
PJx, cr) = Pqgp(o-) [r(x) + 3xr'(x)] , 
Pllix, cr) = Pqgp(o-) [r(x) - 6xr' (x)] . 

where the functions £qgp and Pqgp define the realistic QCD equation of state constructed in 



Ref. 111411 . The function r(x) characterizes properties of the fluid which exhibits the pressure 



r(x) — 



anisotropy x |01 

i-"3 , xarctan 

1 + . (5) 



In the isotropic case x = 1, r(l) = 1, r'(l) - 0, and Eq. ^ is reduced to the standard equation of 
state used in fl4] . 

The function Z in Eq. (|2]l defines the entropy source. We use the form proposed in 10] 

(1 - yi)^ cr 

2(cr,x) = , (6) 

yX Teq 

where the time-scale parameter Teq controls the rate of equilibration. In this work we use the 
constant value Tgq — 1 fm. In Refs. Js, T^, 12] the medium dependent Teq was used, which 



was inversely proportional to the typical transverse momentum scale in the system. If a constant 
value of Teq is used, the system approaches the perfect fluid behavior for t » Tgq. 

In the limit of small anisotropy Eq. © is consistent with the quadratic form of the entropy 
production in the Israel-Stewart theory fSj. Far from equilibrium, hints for the form of S are 



lacking, although we may expect some suggestions from the AdS/CFT correspondence |15|]. 
Thus, for large anisotropics the formula (|6]l should be treated as an assumption defining the 
dynamics of the system. 



3. Initial conditions and freeze-out 

In the general 3-1-1 dimensional case we have to solve Eqs. ^ and (|2]) for cr, x, u^, Uy, and §, 
which depend on t,x± - ix,y), and 77 (t is the proper time and 77 is the space-time rapidity). We 
fix the initial starting time to tq = 0.25 fm. Similarly to other hydrodynamic calculations, we 
assume that there is no initial transverse flow, Ux{to,x±, rf) = Mv(to, x^, 77) = and that the initial 
longitudinal rapidity of the fluid is equal to space-time rapidity, i?(to, x^, /;) - rj. In this text we 
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Figure 1: (Color online) Transverse-momentum dependence of the elliptic flow coefficient vt of + 
calculated for c = 20 - 40% {b = 7.84 fm) at midrapidity and for Teq = 1.0 fm/c; xg = 100 (dashed blue 
lines) and xq = 0.032 (dotted green lines). The results are compared to the PHENIX Collaboration data 
(red dots) y^. 



present the results for two scenarios: i) the initial source is strongly oblate in momentum space, 
x(t[),x±,ti) = 100, and ii) the source is prolate in momentum space, x{to, Xj^,ri) = 0.032; the 
latter value is chosen since r(lOO) - r(0.032). The initial entropy density profile has the form 

0-(){j], = Cr(To, ?7, X^) = figqp [fii p(b, 7], xj] , (7) 

where b is the impact parameter, and p(b, ri, x_i) is the normalized density of sources, p{b, rj, x_i) - 



p{b, 7], x±)/p(0, 0,0), for details see 11 1311 . The quantity ei is the initial energy density at the 
center of the system created in the most central collisions. Its value is fixed by the measured 
multiplicity, separately for two different physical scenarios considered in this paper. We use Si - 
48.8 GeV/fm^ and 80.1 GeV/fm^ for xq = 100 and xq = 0.032, respectively. 

The evolution is determined by the hydrodynamic equations until the entropy density drops 
to (Tf - 1.79 fm"^, which for x - I corresponds to the temperature Tf = 150 MeV. According 
to the single-freeze-out scenario, at this moment the abundances and momenta of particles are 
expected to be fixed. The processes of particle production and decays of unstable resonances are 
described by using THERMINATOR | 14,l1Z]^ which applies the Cooper-Frye formalism to generate 
hadrons on the freeze-out hypersurface. 



4. Results 

The model results describing the pseudorapidity distributions, transverse-momentum spectra, 
the elliptic and directed flow coefficients, and the HBT radii have been obtained with different 
initial anisotropics of pressure; xq - 100 (dashed blue lines) and xq - 0.032 (dotted green lines). 
In all of the considered cases we find good agreement between the model results and the data. 
Moreover, we find that the results obtained with different initial anisotropics are practically the 
same. This is so because we have adjusted the initial energy density separately for two different 
values of xq. A larger (smaller) initial energy density is used for the initially prolate (oblate) 
system. Our results describing the elliptic and directed flow are shown and compared to the 
RHIC data (Au-i-Au collisions at the highest beam energy y^NN - 200 GeV) in Figs. 1 and 2. 
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Figure 2: (Color online) Pseudorapidity dependence of the directed flow of charged particles for the centrality 
bin c = 5 - 40% and r^q = 1.0 fm/c; xq = 100 (dashed blue lines) and xo = 0.032 (dotted green lines). 
The results are compared to the experimental data from STAR (red dots) fl^ and PHOBOS (green squares) 

Hi. 



Our results indicate that the final hadronic observables are not sensitive to the early anisotropy 
of pressure. The flows are built up during the whole time evolution of the system, hence the rela- 
tively short early anisotropic stage does not influence the results. In our opinion, the insensitivity 
of flie hadronic observables helps us to circumvent the early thermalization/isotropization puzzle. 
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